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Abstract 

We investigate the Dirac equation in the semiclassical limit — > 0. A semi- 
classical propagator and a trace formula are derived and are shown to be 
determined by the classical orbits of a relativistic point particle. In addition, 
two phase factors enter, one of which can be calculated from the Thomas 
precession of a classical spin transported along the particle orbits. For the 
second factor we provide an interpretation in terms of dynamical and geomet- 
ric phases. 
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The first one to seek a semiclassical treatment of tlie Dirac equation in tlie spirit of the 
WKB-method appears to be Pauh who gave a solution for a special case. He found 
that in the semiclassical limit the translational motion is independent of the spin degrees 
of freedom. Because of this fact the formalism was criticized by de Broglie with the 
remark that one would expect 'classical objects' like electromagnetic moments to influence 
the trajectories. This controversy was clarified by Rubinow and Keller [^] in a paper that 
seems to have been overlooked by some later authors. Rubinow and Keller pointed out 
that the moments of an electron are proportional to h so that in leading order as — *• 
the influence of spin on the trajectories vanishes. However, in next-to-leading order the 
dynamical equation for Thomas precession is obtained from the Dirac equation. Since 
only the ratio of the magnetic moment and spin enters this equation, it contains no h and 
therefore can be interpreted as describing the dynamics of a classical spin. 

The general set-up for semiclassical quantization in the case of multicomponent wave 
equations was developed by Littlejohn and Flynn ||^ . In a short-wavelength approximation 
they replaced the matrix-valued wave operator by a matrix-valued Hamiltonian function, 
such that its eigenvalues generate Hamiltonian dynamics in phase space. But even if these are 
integrable, an application of EBK quantization was found to be obstructed by the presence 
of additional phases. In a formalism was presented that allows to treat matrix Hamil- 
tonians with no (globally) degenerate eigenvalues. Subsequently, Emmrich and Weinstein 
outlined how to proceed in the degenerate case that, e. g., occurs for the Dirac equa- 
tion, and pointed out the problems of formulating a Bohr-Sommerfeld quantization. They, 
moreover, uncovered the global geometric meaning of the additional phases. A semiclassical 
quantization for special configurations, based on the complex WKB-method, is presented in 

i- 

In this paper we will follow an alternative approach in that we investigate the semiclassi- 
cal time evolution and then set up a trace formula. This procedure avoids (some) difficulties 
that one encounters with semiclassical approximations to eigenspinors and, furthermore, 
is not restricted to classically integrable systems. We basically follow the approach that 
was developed by Gutzwiller for the Schrodinger equation. Hence the basic object to 
be studied is the integral kernel K{x,y,t) of the time evolution operator U{t). Gutzwiller 
represented the kernel by a path integral and evaluated this semiclassically. However, here 
we prefer to use a representation of the kernel in terms of an oscillatory integral. This pro- 
cedure can be made mathematically rigorous as, e. g., explained in for the Schrodinger 
equation. In a second step we pass to the energy domain via Fourier transform, and then 
take the trace over spatial coordinates as well as over spin degrees of freedom. This results in 
a periodic orbit formula for spectral functions. Special attention is paid to the role of spin. 
Our philosophy of a systematic semiclassical expansion in the context of the Dirac equation 
automatically ensures that spin is treated quantum mechanically from the outset, without 
any ad hoc semiclassical approximation. As mentioned above, the semiclassical asymptotics 
introduces an adiabatic decoupling of (slow) translational and (fast) spin degrees of freedom. 
This happens in such a way that to lowest orders in h the expected dynamical equations for 
both kinds of degrees of freedom emerge. In addition, our procedure allows to re-interpret 
the additional phases in terms of dynamical and geometric phases associated with a pre- 
cessing spin. The degree of freedom that is lost upon passing from a quantum mechanical 
description of spin in terms of SU(2)-matrices to a classical description in terms of vectors 
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s G with fixed length |s| can be reconstructed from one of these phases. A detailed 
account of our approach will be presented elsewhere [|T0[] . 

Let us now briefly summarize the calculations and results. We investigate the Dirac 
equation 

ih—^ = HD^ix,t) (1) 
with the (quantum) Hamiltonian 

Hn := a (^-V - -A{x)^ + I3mc^ + e (^{x) (2) 

that acts on a suitable domain in the Hilbert space L^(R'^) ® C'^. The Dirac algebra is 
realized by 

«=f":l and : ) , (3) 



where <t is the vector of Pauli matrices. The time evolution kernel is defined by 

v^(x,t)= f K(x,y,t)vI/o(y)d3|/ (4) 

so that it has to fulfill the Dirac equation for t > with initial condition K{x, y, 0) = 
'Q.4x4:S{x — y). Anticipating the occurrence of solutions of appropriate classical equations of 
motion for positive and negative energies, respectively, we choose the semiclassical ansatz 

K{x, y, t) = / [a+e^^^ + a^e'^'^-] d'^ (5) 



{2nh) 



with phase functions (j)^ = (j)^{x, y, t; ^). The amplitudes are 4 x 4-matrices with semi- 
classical expansions 

oo 

a±(x, y, t; = ^^(-i/i)' (^c, y, t; ^) . (6) 

k=0 

In order to account for the initial condition of the kernel, we have to choose 0^|t=o = {x — y)$ 
and a~^\t=o + (^h\t=o = l4x4- Inserting (H) into (|l]) and comparing like orders in h yields to 
lowest order matrix equations which have solutions with non-zero only if 0^ satisfy the 
Hamilton- Jacobi equations 

i/±(V.0±,a.) + ^ = O (7) 
with the (classical) Hamiltonians 



H^{p, x) = eip{x) ±^c^{p- -^M^)) + ^^c4 . (8) 
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These are the (twofold degenerate) eigenvalues of the matrix- valued symbol cx{p — e/cA) + 
Pmc^ + eip of Hd. Due to (|^) one can separate y in 0=*= according to 0^ = S'^(cc, ^, t) — y$. 
When one applies the method of stationary phase to (|^) as /j. — >• 0, it turns out that at 
stationary points 5** generates a canonical transformation that describes the dynamics of a 
relativistic point particle from y to a; in time t. 

We now turn to the equations that occur in next-to-leading order in h, and which contain 
terms involving both and af. Here we restrict to the index +. An equation for only 
is obtained through a multiplication on the left with the hermitian conjugate of the 
4 X 2-matrix 



Vt:= V{x,y,t;^) 



e + mc 

(TTZ 



(9) 



with e := Vc^tt^ + mc^ and tt = x4> ~ f ^) whose columns are the eigenvectors associated 
with . (We will denote by Wt the corresponding matrix of eigenvectors associated with 
H~ .) We then define a 2 x 2-matrix by 



(10) 



and remark that only this construction with Vt on the left, together with the appropriate 
0^, ensures that the equation to lowest order in h is fulfilled. Moreover, at t = the initial 



condition h 



+ |i=0 



Il2x2 ensures that a. 



Q |j=o = is ^he projector on the i?+-eigenspace. 
A respective remark applies to Oq so that the initial condition for + Oq is fulfilled. An 
obvious interpretation of this ansatz is as follows. Given an initial 4-spinor \E'o at time t = 0, 
Vq projects it onto the if"'"-eigenspace and converts it to a 2-spinor. This is propagated 
to time t and then Vt maps the 2-spinor back to the 4-spinor representation. Using the 



ansatz (10) in the equation of next-to- leading order in h then yields the following transport 
equation for b = b+, 



Vpif+(V,5+,x)V, 



d_ 

di 



-(Mi + iM2)6 



(11) 



with the hermitian 2 x 2-matrices 

3 / 3 

2 



Ml 



+ 



ec 



^ dpjdpk dxkdxj dpjdxj 



1 



2x2 



Ma := aB + 

2e 



ec 



2e(e + mc^ 



-cr(7r X E) . 



(12) 



We need to solve (^) only along the orbits in phase space, in which case the left-hand side 
can be viewed as the total time derivative b along these orbits. To arrive at ([TTD we used 
Coulomb gauge, but this doesn't restrict the result because it only contains the fields E and 
B. The contribution to (0) coming from Mi is well known from the Schrodinger case 0, 
and therefore the ansatz 



'det 



(13) 
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with some 2 x 2-matrix d, proves usefuL From (|TT]) one then obtains the transport equation 

d + iM2d = 0, d\t=o = -SL2X2 , (14) 

for d, which only involves the spin degrees of freedom. Due to the unitarity of the time 
evolution and the initial condition, d has to be an SU(2)-matrix. 

The additional phases discussed in are caused by M2. The second term in M2, 
see (0), can be shown to be V"^ pH^'V x + -§i)V^ and thus is a projection of the natural 
connection on the trivial C^-bundle over phase space onto the if"'"-eigenbundle. According 
to |1TT1J6[], it hence is the Berry term identified in P|. We call this SU(2)-Berry term in order 
to distinguish it from the U(l)-phase originally introduced by Berry [|1^. The first term 
in M2 then is the "no name term" of ^ that has been shown to be related to a Poisson 
curvature in 0. In fact, it measures to what extent the classical time evolution tends to 
leave the if^-eigenspace. In physical terms, the first (curvature) term is the interaction 
of spin and magnetic field, and the second (SU(2)-Berry) term represents the spin-orbit 
coupling. Analogous considerations apply to H~ . 

We are now in a position to state the following semiclassical expression for the time 
evolution kernel. 



K{x,y,t) 



(27rin)3/2 



with 



-7* 
Ixy 



\ 



( ^ fxy 



Ixy 



7xy 



* ^ ixy 



det 



ixy 

dxjdyk 



Ixy ^ T 



{i + om 



(15) 



where 7^^^ labels the classical orbits that connect y and x in time t. is Hamilton's 
principal function, which is the Legendre transform of with respect to ^, and z/^ is the 
Morse index of the corresponding orbit. 

We are still left with the calculation of d. Since d G SU(2), we can use the representation 



d 



with 



I |2 I I |2 

\u\ + \v\ 



1 . 



(16) 



A candidate for a 'classical spin' should be a vector s G with fixed length, which we find 
convenient to choose as |s| = 1. We thus seek a map d y-* s from SU(2) to C M^. To 
achieve this we propose to use the well known Hopf map tth '■ SU(2) — > defined by 



s :- 



2Re{uv) 
2 lm{uv) 
InP — If P 



M, f cr 



(17) 



The last equality reveals that s is also connected to a suitable spin expectation value. From 
(p!7|) and (jM]) it follows that s fulfills the classical equation 



ec 



e(e + mc^ 



-IT X E 



\t=o 




(18) 
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describing a precessing spin. After [Q, this is commonly called Thomas precession, see also 
T3| . Due to the initial condition, s G will stay on the northern hemisphere for sufficiently 



small times. We then choose polar coordinates. 



sin^^ cos( 

sin 6' sin0 I , (19) 



cos 6* 

which allows to calculate c? up to a phase 77, where u/\u\ = e'^. 



, cos(^^/2)e"' -sin(^/2)e-(^-<^) ^ . 
^ sin(0/2)e*(''-'^) cos(^/2)e-*'' ' ' ^ ^ 



The equation for rj, which is obtained upon inserting (pO]) into (|T^ and multiplying by d\ 
can immediately be integrated, 

V=l f s(-B--^^7rxE]dt' + l Al-cos0)0dt'. (21) 

The first term is a dynamical phase associated with the energy of a (classical) magnetic 
moment in given electromagnetic fields, whereas the second term is a geometric phase. We 
remark that once s enters the southern hemisphere of S^, one should change the phase 
convention in that v/\v\ = e'^^ is used to describe the non-classical degree of freedom. In 
(PI) this amounts to replace 1/2(1 — cos^^)d0 by —1/2(1 + cos6')d0. These two expressions 
are the well known gauges of the vector potential for a magnetic monopole of strength 1/2 
situated at the origin of the sphere. The geometric phase caused by this connection is 
reminiscent of (but not identical to) the quantum mechanical Berry phase of a precessing 
spin [|T2 |. 



Now all terms appearing in the semiclassical time evolution kernel are fixed. A non- 
relativistic approximation is obtained, if in (^) one only keeps the leading asymptotic term 
as c — > 00. As a result, one is left with a block diagonal formula. On the other hand, we 
also performed the above program of a systematic semiclassical expansion in the case of the 
Pauli equation. Its result coincides with the upper left block of the former approximation. 

Our next goal is to derive a semiclassical trace formula from (|15|). Since always has a 
continuous spectrum, which contains at least (—00, —mc^) U (mc^, 00), we find it convenient 
to introduce an energy localization such that finally only the discrete spectrum of Hd enters. 
We thus assume that the spectrum of Ho is purely discrete on an interval I = {Ea, E^). Then 
we choose a smooth function x{E) which is non-zero only on J, such that x(-^n) = 1 for all 
eigenvalues En- This can always be achieved if there is no accumulation of eigenvalues at Ea 
or Ely. Instead of the full time evolution operator we then study its restriction x{HD)U{t). 
To leading order in this restriction only causes additional factors xi^-y^y) in (plSj). For 
Schrodinger operators this procedure is described in The restricted time evolution kernel 
has a spectral representation 

K{x, y,t) = J2 xiEn) '^n{x)<fi{y) e-i^"* (22) 

n 

with orthonormal eigenspinors We define a regularized Green function by 
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G<^{x,y,E):-- 



1 

2^ 



g{t)e^^^ K{x,y,t) dt 



(23) 



where ^ is a smooth test function such that its Fourier transform g vanishes outside a finite 
interval. Taking the trace of over spatial variables and matrix components yields 



Tr G^) (E) = Tr4x4 £^ G^^ix, x, E) d?x = x{En) Q (^%-^) • 



(24) 



The trace formula can now be derived from (23) and ( plf ) when one introduces the semiclas- 
sical approximation (^), but now modified as described above in order to apply to the kernel 
K. As in the case of the Schrodinger equation, the integrals necessary to calculate (0) can 
be evaluated with the method of stationary phase. The first and foremost contribution then 
derives from the stationary points with t = 0. In leading semiclassical order this term (also 
called Weyl term) involves the volumes \^%\ of the energy shells in phase space, 



6{H^{p,x) - E) d^pd^x 



(25) 



Up to terms 0{h°°), all further contributions are caused by the non-trivial periodic orbits 
of the classical dynamics generated by if+ and H^. In case that all periodic orbits are 
isolated and unstable (i. e., hyperbolic or inverse hyperbolic) these contributions will be 
given explicitly. Their calculation is exactly parallel to the case of the Schrodinger equation. 
The only additional factor that enters comes from the trace over the spin degrees of freedom. 
If T is the period of a periodic orbit, Vp = Vq so that 

Tr2 X 2 (K|VT(i+ 



Tr4x4(V^Trf+K)^) 

We now choose E E I such that x{E) 



- Tr2x2'^+ = 2cos(6'/2) cos?7 . 
1, and thus obtain the trace formula 



(26) 



h 



g(o) + \n} 



27r 



{l + 0{h)} 



7? 



2n 



A^^e"''-^^''^-''^^'-t {l + 0{h)} 



with A + 

IP 



2Tt cos( 

7p 



^^±/2)cosr/^± 



det(M^± 



(27) 



On the right-hand side the sum extends over the classical periodic orbits 7^ of energy E. 
Furthermore, S{E) = ^ pdx is the action, T the period, /i the Maslov index, and M is the 
(linearized) Poincare map; T* denotes the associated primitive period. 

The factor 2 cos(6'/2) cos// emerging from the spin degrees of freedom has to be inter- 
preted as follows. The angle 6 measures the discrepancy between the directions of the spin 
vector after this has been transported along a given periodic orbit with the dynamics dic- 
tated by (0). The contribution of the periodic orbit to the trace formula is then weighted 
with cos(0/2). The second term arises from quantum mechanics and, as explained above, is 
composed of a dynamical as well as of a geometric phase. The factor of two finally indicates 
the presence of two spin directions. 
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